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Abstract 

There has been overwhelming evidence that coherent structures play 
a critical role in determining the overall transport in a variety of sys- 
tems. We compute the probability distribution function (PDF) tails 
of momentum flux and heat flux in ion-temperature-gradient turbu- 
lence, by taking into account the interaction among modons, which 
are assumed to be coherent structures responsible for bursty and in- 
termittent events, contributing to the PDF tails. The tail of PDF 
of momentum flux R = (v x v y ) is shown to be exponential with the 
form exp {— £i? 3 / 2 }, which is broader than a Gaussian, similarly to 
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what was found in the previous local studies. An analogous expres- 
sion with the same functional dependence is found for the PDF tails 
of heat flux. Furthermore, we present a detailed numerical study of 
the dependence of the PDF tail on the temperature and density scale 
lengths and other physical parameters through the coefficient £. 
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I Introduction 



One of the main challenges in magnetic fusion research has been to predict 
the turbulent heat and particle transport originating from various micro- 
instabilities. The ion-temperature-gradient (ITG) mode is one of the main 
candidates for causing the anomalous heat transport in core plasmas of toka- 
maks pQ. Significant heat transport might however be mediated by coherent 
structures such as streamers and blobs through the formation of avalanche 
like events of large amplitude [2]- [I], as indicated by recent numerical stud- 
ies. These events cause the deviation of the probability distribution functions 
(PDFs) from a Gaussian profile on which the traditional mean field theory 
(such as transport coefficients) is based. In particular, PDF tails due to 
rare events of large amplitude are often found to be substantially different 
from Gaussian although PDF centers tend to be Gaussian |5|. These non- 
Gaussian PDF tails are manifestations of intermittency, caused by bursts and 
coherent structures. The characterization of these PDF tails thus requires a 
non-perturbative method. 

There are many coherent structures of interest, including zonal flows and 
streamers. While streamer-like structures (ke 3> k r ) enhance transport zonal 
flows can dramatically reduce transport [6]- [Tj. The zonal flows are poloidally 
and toroidally symmetric (kg = 0, kn = 0) and radially inhomogeneous (k r ^ 
0) flow structures in toroidal plasmas. It would thus be of great importance 
to develop a theory of the formation of coherent structures and PDFs of heat 
flux due to these structures. 

The purpose of this paper is to investigate the likelihood of the forma- 
tion of coherent structures by computing the PDF (tails) of the Reynolds 
stress and predict the PDFs of heat transport. Specifically, we extend a non- 
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perturbative theory of the PDFs of local momentum flux and heat flux to ac- 
commodate global fluxes by incorporating the interactions among structures. 
An advanced fluid model for the ITG mode is used [8] that has been successful 
in reproducing both experimental [H] and non-linear gyro-kinetic results [ID] . 
The advanced fluid model is expected to give qualitatively and quantitatively 
accurate results in tokamak core plasmas or the flat density regime. For in- 
stance, this was shown in the plot of the slope of the ITG threshold as a 
function of e n in the flat density regime (see Figure 1 in Ref. (TTJ) where this 
model recovers the full linear kinetic result within 5%, far better than the 
prediction from an approximate kinetic model using a constant energy ap- 
proximation in the magnetic drift. In particular, the gyrokinetic simulations 
of ITG mode turbulence reported in the Cyclone work Ref. [10] indicate that 
there is a strong excitation of zonal flows close to marginal stability where the 
non-linearly generated flows were able to damp out the turbulence resulting 
in a non-linear up-shift in the critical temperature gradient needed to obtain 
transport for longer time scales. One particular coherent structure in ITG 
drift wave turbulence is the two dimensional bipolar vortex soliton solution 
called the modon [12]- [2]. The modon is an exact solution to the non-linear 
governing equations for ITG mode turbulence and travels perpendicular to 
both the strong magnetic field and the background density gradient. The 
PDF tail is viewed as the transition amplitude from an initial state with no 
fluid motion to a final state governed by the modon with different amplitudes 
in the long time limit. 

The theoretical technique used here is the so-called instanton method, a 
non-perturbative way of calculating the PDF tails. The PDF tail is first for- 
mally expressed in terms of a path integral by utilizing the Gaussian statistics 
of the forcing. An optimum path will then be associated with the creation 
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of a modon (among all possible paths) and the action is evaluated using 
the saddle-point method on the effective action. The saddle-point solution 
of the dynamical variable <p(x,t) of the form 4>(x,t) = F(t)ip(x) is called 
an instanton if F(t) = at the initial time and F(t) ^ in the long time 
limit. The instanton is localized in time, existing during the formation of 
the modon. Thus, the bursty event can be associated with the creation of a 
modon. Note that, the function ip(x) here represents the spatial form of the 
coherent structure. Historically, the instanton method was used in gauge field 
theory for calculating the transition amplitude from one vacuum to another 
vacuum [15]- [16] . 

We will show exponential PDF tails for both heat flux (H) and momentum 
fluxes (R) of the form e~^ H3/2 and e~^ R3/2 , similar to the PDFs of local fluxes 
in the earlier works [17]- [19]. The dependence of the coefficient £ (i.e. the 
overall amplitude of the PDF tails) on parameter values (e.g. the density 
gradient, temperature gradient, curvature, etc) will be studied in detail. We 
will also demonstrate the generation of monopole vortex via the interaction 
of two dipoles (modons). 

The paper is organized as follows. In Sec. II the model for the ITG mode 
is presented together with preliminaries for the path-integral formulation for 
the PDF tail of momentum flux. In Sec. Ill the instanton solutions are 
calculated and the PDF tail of momentum flux is estimated in Sec IV. We 
provide numerical results in Sec. V and a discussion of the results and the 
PDF tail of heat flux and conclusion in Sec. VI. 
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II Non-perturbative calculation of momentum 
flux PDF 



The ITG mode turbulence is modeled using the continuity and temperature 
equation for the ions and considering the electrons to be Boltzmann dis- 
tributed; quasi- neutrality is used to close the system [8] . In the present work 
the effects of parallel ion motion, magnetic shear, trapped particles and finite 
beta on the ITG modes are neglected since in previous works that the effect 
of parallel ion motion on the ITG mode was found to be rather weak [20] . 
The continuity and temperature equations are, 
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Here [A, B] = (dA/dx)(dB / dy) — (dA/ dy)(dB / dx) is the Poisson bracket; / 
is a forcing; n = (L n J p s )5n/n , <j> = (L n J p s )e5(p/T e , f 4 = (L n / p s )5Ti/T i0 are 
the normalized ion particle density, the electrostatic potential and the ion 
temperature, respectively In equations (1) and (2), r = Tj/T e , p s = c s /Q C i 
where c s = ^jT e /rrii, Q C i = eB/rriiC and v is collisionality. We also define 
L f = - (dlnf/dr)" 1 (/ = {n,^}), m = L n /L Ti , e n = 2L n /R where R is 
the major radius and a« = r (1 + 7)i). The perpendicular length scale and 
time are normalized by p s and L n /c s , respectively. The geometrical quanti- 
ties are calculated in the strong ballooning limit (8 = 0, gi (8 = 0, k) = 1/ k 



where gi (8) is defined by ujd (ff) = u)*e n gi {8)) [21]- [22], with a;^ 



p s c s k y /L n . The system is closed by using quasi-neutrality with Boltzmann 
distributed electrons. Note that the approximate linear solutions to Equa- 



tions (1) and (2) with / = give the dispersion relation with real frequency 
(u r ) and growth rate (7) as, 

u r = n jz ky , 2 r (l ~ ( 1 + ^r) e n9i ~ fcl (a» + \re n gX] , (3) 
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Finite Larmor Radius (FLR) effects on the rjuh are here neglected, although 
these are included in the numerical study leading to Figures 1 and 2. 

We formally calculate the PDF tails of momentum flux by using the 
instanton method. To this end, the PDF tail is expressed in terms of a 
path integral by utilizing the Gaussian statistics of the forcing / [23]. The 
probability distribution function for Reynolds stress R can be defined as 

P(R) = (6((v x v y ) - R)) 

dX exp(iXR) (exp(—iX(v x v y ))) 

dXexp(iXR)I\, (6) 



where 

h = (exp(-iXv x v y )) . (7) 
The integrand can then be rewritten in the form of a path-integral as 

V4>V(j)e' s \ (8) 



The angular brackets denote the average over the statistics of the forcing /. 
By using the ansatz Tj = \4>i which will be justified later (see Eq. (15)), the 
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effective action S\ in Eq. (8) can be expressed as 

Sx = -./^0^-(^-« i ^)Vi0+(l-6„^)g-/3[0,Vi0]) 
+ -J dtd 2 xd 2 x'(j){x)K(x — x')4>{x') 



Here, 



(3 = 1 + r + rx, (10) 
rji - |(1 - C/) 

* = 77 J L n • 11 

In Eq. (11), U is the modon speed (see Eq. (15)). To obtain Eq. (9) 
we have assumed the statistics of the forcing / to be Gaussian with a short 
correlation time modeled by the delta function as 

(f(x,t)f(x',t')) = 5(t-t')K(x-x'), (12) 

and (/) = 0. The delta correlation in time were chosen for the simplicity of 
the analysis. In the case of a finite correlation time the non-local integral 
equations in time are needed. We will also make use of the completeness 
of the Bessel function expansion and write k(x — x') = K (J (kx) Jo(kx') + 
J l (kx)J 1 (kx')(cos 9 cos tf'+sin 9 sin 6')+J 2 (kx)J 2 (kx') (cos 29 cos 20'+sin 29 sin 29') + 
...). 



Ill Instanton (saddle-point) solutions 

We have now reformulated the problem of calculating the PDF to a path- 
integral as in Eq. (6). Although the path integral cannot in general be 
calculated exactly, an approximate value can be found in the limit A — > oo by 
using a saddle point method. The idea of the saddle-point method is that the 
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integrand has a unique global maximum and that all significant contributions 
to the integral come only from points in the vicinity of this maximum. In 
the limit A — > oo a particular path that satisfies the saddle-point equations 
gives the leading order contribution. The saddle-point equations are; 

t = °' (13) 
f - ° <"> 

Since a direct application of the saddle-point equations results in very compli- 
cated partial differential equations for and 0, we assume that the instanton 
saddle-point solution is a temporally localized modon. That is, we assume 
that a non-linear vortex soliton solution exists to the system of Eqs (l)-(2) 
by assuming that the electric potential can be written 

<i>{x,V,t) = ^(x,y- Ut)F(t), &nd Ti = x<fi- (15) 

The function ip in Eqs (16)-(17) is the spatial form of the coherent structure, 
which is assumed to be the sum of the two modons with the ratio of strength 
e as follows, 

oi 

ip(x,y — Ut) = ciJi(kr)(cos9 + e sin 60 + — r cos 9 for r < a, (16) 

k 

ip(x,y — Ut) = C2-Ki(pr)(cos6 l + e(r) sin 9) for r > a. (17) 
Here J\ and K\ are the first Bessel function and the second modified Bessel 



function, respectively. Here r = y/x 2 + y 2 , tan 9 = y'/x, y' = y — Ut, a = 
(At - k 2 A 2 ), A x = (1 - e n gi - U)/f3, A 2 = (U + oa)/f3. By matching the 
inner and outer solution at r = a we find the conditions c\ = —aa/Ji(ka), 
c 2 = -Ua/Kxipa), J[(ka) / 'J x {ka) = (1 + k 2 /p 2 )/ka - kK[(pa)/pKi(pa); 
U is the velocity of the modon, and a is the size of the core region. The 
function e(r) is chosen such that the matching conditions are similar to those 
in previous previous studies [17]- [18]. 



It is important to note that if e = that the averaged Reynolds stress 
vanishes. In general the Reynolds stress associated with the modon can be 
found to be 

d / v f A 2 , dtp dtp. 

= -errc? y drr[— J (fcr) + ^(^r^) + 2 K ) 7 

+ ^ViW^ + ^^CAt)) 8 ]. (18) 

The Reynolds stress in Eq. (18) represents the effective force driving coher- 
ent structures. Interestingly, it shows that the non-linear interaction of two 
modons (dipoles) can generate a monopole, given by the zeroth order Bessel 
function (Jo), as well as other more complicated structures. Given that a 
monopole is the dominant term our result, Eq. (18), indicates the tendency 
of monopole formation from modons. This tendency was observed in numer- 
ical studies in Ref. [H], [21]. Next, the action S\ is to be expressed only as 
an integral in time by using the conjugate variables 



F = J d 2 x0(x,t)J o (kr), (19) 

F u = J d 2 x<p{x,t)J 1 (kr)sm6, (20) 

F lc = J d 2 x(j)(x,t)J 1 (kr)cos6, (21) 

F 2s = J d 2 x0(x,t)J 2 (kr) sin2#, (22) 

F 2c = f d 2 x(p{x,t)J 2 (kr) cos 26. (23) 



Note that the contribution from the outer solution (r > a) to S\ is neglected 
compared to that from the inner solution (r < a) for simplicity. The outer 
solution decays fast and inherently gives a minor contribution to the PDF 
tail. The action S\ consists of three different parts; the ITG model, the 
forcing and the Reynolds-stress parts respectively. The ITG model part of 
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the action can be reduced to 
Iitg — — i J d 2 xdt 



F((l + k 2 )ip — k 2 ax) 
+ F[l-e ngi p-U-k 2 (U + ai)+pFk 2 a 
+ Fuk 4 ^ — Fuk 4 ax 



chb 
dy 



(24) 



The full action including the forcing and Reynolds stress terms can then 
be expressed in terms of F, F and the conjugate variables F 

S x = -ij dt[nF(F lc + eF u ) + F{ l2 F 2s + e 73 F + e 74 F 2c ) 

+ F 2 { lb F 2s + e 76 F + e l7 F 2c ) + 7s (F lc + eF u ))\ 

+ 1 -k J dt(F 2 + 2(F 2 C + F 2 S ) + 2(F 2 S + Fl)) 

- XR J dtF 2 5(t) (25) 

Here the coefficients are; 

(26) 
(27) 
(28) 
(29) 
(30) 

H). (31) 

The equations of motions for the instanton are found by the variations of 
the action with respect to F, F , Fi c , Fi s , F 2c and F 2s ; 
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(32) 
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-te{ l3 F + 7e F 2 ) + k F = 0, (33) 
= -i( 7l F + 78 F) + = 0, (34) 

= -ie(F + 7s F) + 2k F 1s = 0, (35) 



s 



-*e( 74 F + 1 7 F 2 ) + 2k F 2c = 0, (36) 



SF 2c 

-K12F + 75 F 2 ) + 2k q F 2s = 0. (37) 



<5F 2s 

The equation of motion for F is derived for t < using Eqs. (32)- (37) as; 
1 <iF 2 

- 7l 2 (l + e 2 )^ = Vl F + 3 V2 F 2 + 2 V3 F 3 , (38) 

77! = 7 2 + 2e 2 7 3 2 + e 2 7 2 + 7 2 + e 2 7 2 , (39) 

V2 = 7 2 7 5 + 2e 2 7 3 76 + e 2 7 4 77 , (40) 

V3 = 75 2 + 2e 2 7 2 + e 2 7 2 . (41) 

The contribution from the dissipation {y) to the term involving the time 
derivative of F (F) cancels out and the equation of motion is exactly solvable. 
In the limit of A — > 00 the relation can be written, 



F = \Z^rbf 2 (42) 

The initial condition is found by integrating Eq. (32) over the interval [—5, 0] 
((5 < 1) and observing that the conjugate variables mediating between the 
forcing and F vanish for t > 0. This can be interpreted as a "causality" 
condition 

i7i(Ac(-<5) + eF ls (-5)) + 2Ai? F(0) = 0. (43) 

The elimination of F lc and F ls using Eqs. (34) and (35) gives a relationship 
between F and F in the limit A — > 00 as, 

7l 2 (l + t 2 )F(-5) ~ 4k \R F(0). (44) 
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Finally, Eqs. (43) and (44) gives the initial condition for F 

F(0) - . (45) 

VW(l + e 2 ) 

We use this initial condition to compute the saddle-point action and then 
predict the scaling of S\ (as A — > oo) to compute the PDF tail in the next 
section. 



IV The PDF tail 

The PDF tail is found by calculating the value of S\ at the saddle-point. Of 
particular interest is the dependency of the action on A in the limit A — > oo 

which will give the PDF tail; 

/ ^ 2 

1 / AXkqRq 

S x ^ --Ai? ; = 

3 VvW(i + e 2 ) 

~ -^h\ 3 . (46) 

The PDF tail of the Reynolds stress (R) can now be found by performing 
the integration over A in Eq. (6) using the saddle-point method, 

P(R) ~ exp(-e(^) 3/2 ), (47) 
2 \/wyKl + e 2 ) 

« - I l tt - ' 48 » 
Here, it is important to note that our assumption that A — > oo corresponds 
to -R — > oo. Equation (47) gives the probability of a Reynolds stress R, 
normalized by the Reynolds stress Rq due to the the modon solution Rq 
given by Eq. (18) which is fixed for given parameters. We have assumed 
that the modon is created and that F(t) = as t — > — oo, which means we 
can interpret the Eq. (47) as a transition amplitude from an initial state, 
with no fluid motion, to final states with different values of R/ Rq. It is also 
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important to note that even though we have assumed a Gaussian forcing 
the tail is non-Gaussian exhibiting intermittency. The overall coefficient £ 
in the PDF in Eqs. (47)- (48) depends on several physical parameter values. 
First, £ — > oo (i.e. PDF vanishes) as the the forcing disappears (k — > 0); 
the instanton cannot form and the PDF vanishes (P(R) — > 0). Second, 
£ increases for larger e leading to a reduction in the PDF. Furthermore, £ 
depends on ion temperature gradient (rji), density gradient (e„), temperature 
ratio (r = Tj/T e ), modon size (a), modon speed (U) and wave number (k). 
These dependencies will be studied in detail numerically in Sec. V. 

V Results 

We have presented a calculation of the tail of the Reynolds stress probability 
distribution function (PDF) in forced ITG turbulence. We have shown by a 
non-perturbative calculation (instanton calculus) that a coherent structure 
can lead to intermittent vorticity flux. A path-integral formulation is devel- 
oped for the PDF tail. The system is solved by assuming that the tail is 
associated with bursty events or the creation of modons (a bipolar vortex). 
The integrals are then estimated by the saddle-point method, which gives 
the functional dependence of the tail as exp {— £(R/ R ) 3 / 2 }. This exponen- 
tial form seems to be ubiquitous in drift wave turbulence. In this section, 
the parameter dependencies of £ will be studied in detail and compared to a 
normalized Gaussian profile. 

In Figure 1, the PDF tails are shown as a function of Reynolds stress 
for forced ITG mode turbulence (blue line), forced Hasegawa-Mima (HM) 
turbulence (red line), Gaussian distributions with same parameters as for 
the ITG mode turbulence (green line) and cases with negative modon speed 
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in ITG mode turbulence (black line) and HM turbulence (dashed black line). 
The other parameters are rji = 4.0, r = 0.5, e n = 1.0, $ = 1, a — 2, U = 2.0 
(U = —5.0 black line and dashed black line), k = 3.0, e = 0.1 and k ps 1.84 
(ITG case), k = 0.81 (ITG with reversed modon speed), k « 1.73 (HM case) 
and k = 1.56 (HM with reversed modon speed). The PDF tail in HM model 
is obtained by setting the parameters (3 = 1 and e n = 0.0. This is equivalent 
to letting r]i = 0.0 and r = Ti/T e = 0.0. It is shown that for U > the PDF 
tail for ITG turbulence is significantly lower compared to the value found for 
HM turbulence, although they are qualitatively similar. The reason for this 
is that the model dependent factor £ is significantly different in ITG mode 
and HM turbulence. In the case of reversed modon speed (U < 0), x changes 
sign and significantly enhances the PDF tail. The enhanced ITG PDF tail 
for the reversed modon speed is due to the change in k; a smaller value of k 
combined with a large negative modon speed gives the enhancement. 

Figure 2 shows the PDF tail as a function of Reynolds stress in ITG 
mode turbulence by varying e n parameter values; k ps 1.91, e n = 0.1 (black), 
e n = 1.0, (blue line), e n = 2.0, (red line) and a Gaussian distribution with 
the same £ as for the e n = 1.0 case (green line). The parameters are U = 4.0 
with all the others the same as those in Figure 1. The PDF tails for different 
density profiles are qualitatively similar, however for peaked density profiles 
(small e n ) the PDF tails are significantly enhanced. 

In most cases, the PDF tails of drift wave turbulence differ significantly 
both quantitatively and qualitatively with the Gaussian distributions. 
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Figure 1: (Color online). The PDF tails are shown as a function of Reynolds 
stress for ITG mode turbulence (blue line), forced Hasegawa-Mima (HM) 
turbulence (red line), Gaussian distributions with same parameters as for 
the ITG mode turbulence (green line) and cases with negative modon speed 
in ITG mode turbulence (black line) and in HM turbulence (dashed black 
line). The other parameters are rji = 4.0, r = 0.5, e n = 1.0, $ = 1, a — 2, 
U = 2.0 (U = —5.0 black line and dashed black line), k = 3.0, e = 0.1 and 
k w 1.84 (ITG case), k = 0.81 (ITG with reversed modon speed), k w 1.73 
(HM case) and k = 1.56 (HM with reversed modon speed). 
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Figure 2: (Color online). Shows the PDF tail as a function of Reynolds 
stress in ITG mode turbulence by varying e n parameter values; e n = 0.1, 
k ~ 1.91 (black), e n = 1.0, (blue line), e n = 2.0, (red line) and a Gaussian 
control distribution with the same £ as for the e n = 1.0 case (green line). 
The parameters are U = 4.0 and the others as in Figure 1. 
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VI Discussion and conclusions 



The present calculation of the PDF tail of momentum suggests the PDF 
tail ~ e~^ R3/2 to be ubiquitous in drift wave turbulence. The 3/2 exponent 
in Eq. (47) follows from the A — ■> oo dependence of the action as S\ ~ 
A 3 and inherently comes from the quadratic non-linearity in the dynamical 
equations. This scaling can be found by balancing the terms in S\ as A — > oo 
as Xcjxp ~ (fxp ~ T<f)<f) 2 ~ from which it follows: <f) ~ A, ~ A 2 , T ~ 

A -1 , and Sa ~ A 3 . Here, T is the typical timescale of the instanton. The 
coefficient £ in Eq. (47) contains all the model dependent information and 
gives the difference between the ITG and the HM models; we have shown that 
PDF tails in drift wave turbulence significantly deviate from the Gaussian 
distributions in most cases. Moreover, it was found that reversed modon 
speed may have significant influence on the PDF tail (e.g see Figure 1). There 
is some interesting recent experimental work done at CSDX at UCSD, where 
the same qualitative Reynolds stress PDF scaling was found (e~^ 3/2 ) [27], 
which agrees with our prediction. 

We now consider the PDF of ion heat flux (TiV x ) given by, 

P(H) ~ (6((T lVx ) - H)} 

= J d\exp{iXH)(exp(-iX(TiV x ))), (49) 

By using the same methodology, we can show that P{H) takes exactly the 
same form as for P(R); i.e. P(H) ~ e~^ H ^ H °^ 3/2 . Here, Hq = (Ttv x ), is 
the heat flux associated with modons. Therefore, heat transport can also be 
significantly enhanced over the Gaussian prediction due to modons. 

In using the instanton method we had to assume the spatial form of the 
coherent structure, a modon in this case, to be an exact solution to the 
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non-linear dynamical equation with fixed parameter values. This means that 
the mechanism of formation of the structure itself was not addressed in this 
analysis. The computation of structure formation will be addressed in a 
future publication. 

We note that for simplicity, we have not included the effect of zonal flows 
in the present study. The system of Eq. (l)-(2) describes the non-linear evo- 
lution of the ITG mode in the presence of a Gaussian forcing (/) , by assuming 
the Boltzmann response. The influence of zonal flows on the PDFs can be 
calculated by using a similar method, and is expected to give a qualitatively 
similar result (e~^ Ri/2 ), but with a different value of the coefficient £. As long 
as the main non-linearity is quadratic in Eq. (1), the same power law should 
follow since the power-law is determined by the highest nonlinearity in the 
equations. Due to the neglect of zonal flows the result obtained in this pa- 
per would be more relevant to understanding transport in drift-wave systems 
that have weak influence of zonal flows e.g. electron-temperature-gradient 
(ETG) turbulence. However, even if zonal flows play an important role in 
the regulation of ITG turbulence, it is still of great interest to investigate the 
PDF tails in ITG turbulence which is weakly regulated by zonal flows. It 
is because many transport simulations have shown that the zonal flows have 
strong effect close to the critical gradient whereas far away from this critical 
gradient the effect of zonal flows is much more complicated. 

In general, for calculating the PDF tail a weighted sum over various 
coherent structures is needed. At present, the only known exact solution is 
the modon which we have assumed to be the underlying coherent structure. 
A generalization should be straight forward if more non-linear solutions were 
available. 

We note that a non-Gaussian scaling of the PDF (the exponent of R) is 
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found even when the forcing is Gaussian, although the exact exponent may 
depend on the temporal and possibly spatial correlation of the forcing (/). 
In the present paper, the forcing is chosen to be temporally delta correlated 
for simplicity. The exponent may also change, if another spatial coherent 
structure is introduced i.e. another non-linear solution to the Eq. (l)-(2) is 
found. The case where the forcing is non-Gaussian (/) will be addressed in 
a future publication. 

Although there are very few numerical simulations of event-size distribu- 
tions, there are however some recent gyrokinetic numerical simulations of the 
spectrum properties of heat flux [25]- [26]. However, the prediction of the 
full spectrum behavior of heat flux is out of the scope of the present work 
and will be addressed in a future publication. 

In summary, this paper presents the first calculation of PDF tails of 
momentum flux and heat flux, which were shown to be significantly enhanced 
over the Gaussian prediction. This suggests that considerable transport is 
mediated by rare events of high amplitude. This is because even if the 
PDF tails have a low amplitude (rare events), these rare events are of high 
amplitude, possibly carrying significant transport. Since the PDF tails are 
enhanced over the Gaussian prediction in our case, these events are more 
likely to mediate considerable transport. The main point of the present 
work is to investigate the influence of different parameters on the coefficient 
(£). The results should be interpreted that for certain parameters, large scale 
events are more likely to be the main cause for transport while less important 
for other parameters. In all cases with parameters relevant for a tokamak 
plasma, the (enhanced) non-Gaussian PDF tail is of great importance. The 
overall amplitude is shown to be larger in ITG than in HM turbulence for 
reversed modon speed (U < 0). 
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